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C"^ . Abstract 

m 

CN ' Starting from a path integral representation of appropriate 4-point and 

^^ ■ 2-point gauge invariant Green functions and from the "Modified Area Law" 

Q>^ ! model , a qq Bethe-Salpeter like equation and a related Schwinger-Dyson equa- 

0^ \ tion can be obtained. From such equations an effective relativistic Hamilto- 

nian can be derived by standard methods and then applied to the determina- 

Qh, tion of the meson spectrum . The entire known heavy-heavy and heavy-light 

Q^l spectra and the lowest light-light Regge trajectories are rather well reproduced 

^ ' in terms of four parameters alone, the light quark masses being fixed a priori 

on typical current values. 

K^ ■ 

X 

^ : 1 Introduction 

In this paper I want to revue a Bethe Salpeter formalism in QCD, which has been 
have developed recently in Milano. Such formalism, if not completely derived from 
first principles, rests, however, only on some non controversial assumptions on the 
Wilson loop correlator |]l|. It generalizes a method introduced previously for the 
case of the heavy quark potential f^ and takes advantage of appropriate Feynmann- 
Schwinger like path integral representations for the QCD Green functions. Its 
three dimensional reduction has been recently applied to the light-light and heavy- 
light quark-antiquark spectrum with very encouraging results Q]. Beside the heavy 
quarkonia, the lowest Regge trajectories for the triplet mm, ms, ss states are very 
well reproduced in slope and intercepts and the known spin averaged light-heavy 
spectrum is obtained up to a mean deviation of about 10 MeV. Notice that it is 





very difficult to obtain a similar result in the frame of a single model. In fact in 
our case the strong coupling constant as, the string tension a and the heavy quark 
masses m^ and rric were already completely determined by the potential fits (apart 
the possibility of a very small rearrangement), while the light quark masses have 
been fixed a priori on typical current values (m^ = rrid = lOMeV, rris = 200 MeV). 
In particular the heavy-light sector is completely parameter free. 

The basic objects from which we start are the ordinary gauge invariant 4-point 
and 2-point Green functions G^\xi,X2,yi,y2) and G^\x — y), to which we refer as 
the "first order" functions. To these ffist order functions certain "second order" ones 
H^\xi, X2,yi, 1/2) and H^\x — y) can be related. It is for the second order functions 
that the mentioned path integral representations can be constructed. 

The important aspect of the above representations is that the gauge field occurs 
simply trough Wilson correlators 



W\r] = -{TiP exp{ig i dx^A^}) . 



associated to loops F made by quark or antiquark world lines and "Schwinger 
strings". In principle these correlators should determine the whole dynamics. Un- 
fortunately, due to confinement and the consequent failure of a purely perturbative 
approach, a consistent analytic evaluation of W from the Lagrangian alone is not 
possible today. However combining incomplete theoretical arguments and lattice 
simulation information various reasonable models can be attempted. 

The most naive but at the same time less arbitrary assumption consists in writing 
i In VT as the sum of its perturbative expression and an area term (modified area 
law (MAL) model) 

i\nW = i{\nW)pert + cTS^in, (2) 

where the ffist quantity is supposed to give correctly the short range limit the second 
the long range one. 

Eq. (!]) is our starting point. In principle any more sophisticated model could be 
used in the context, at the condition that it preserves certain general properties of 
functional derivability of the definition (|I]). In practice not even @ can be treated 
exactly. Actually we shall replace the minimal surface S'min by what can be called 
its "equal time straight line approximation" , which we shall explain later. 

In connection with (Q) it is convenient to consider a third type of Green functions 
H{xi, X2, yi, 2/2) and H{x—y) which are obtained from H^\xi, X2, yi, 2/2) and H^\x — 
y) by omitting in their path integral representation contributions related to the 
Schwinger strings. In the limit of vanishing xi — X2, yi — y2 or x — y such new 
quantities coincide with the original ones and are completely equivalent for what 
concerns the determination of bound states, condensates, chiral symmetry breaking, 
etc. 

For H{xi,X2,yi,y2) and H{x — y), an inhomogeneous Bethe-Salpeter equation 
and a Dyson-Schwinger equation, respectively, can be derived in the configuration 
space, with kernel obtained as an expansion in ctg and a. Such equations can be 



rewritten in a more conventional form in the momentum space and as such apphed 
e. g. to the spectrum of the mesons. To this aim in principle one should solve 
the DS equation first and use the resulting propagator in the BS equation. In 
practice a direct treatment of their full four dimensional expressions seems to be 
a very difficult task. However, if one neglects the pseudo-scalar mesons, a three 
dimensional reduction of the BS-equation seems to be appropriate. This can be 
obtained by standard methods in the form of an eigenvalue equation for an effective 
squared mass operator or simply a relativistic Hamiltonian |]l], H . 

For the case of the pseudo-scalar mesons a complete four- dimensional treatment 
would be imperative. In fact the use of the free quark propagator implied in the 
three dimensional reduction can not be suitable in this case, for the strict interplay 
existing between zero mass Bethe-Salpeter wave function and the propagator in the 
chiral limit; we shall limit to few comments. 

It should be mentioned that our formalism is strictly related from different point 
of view to to the works of ref. |^. 

The plan of the paper is the following one. We discuss the gauge invariant Green 
functions and their path integral representations in sect. 2; modified area law model 
and straight line approximation in sect. 3; Bethe-Salpeter, Dyson-Schwinger equa- 
tions and the problem of chiral symmetry in sect. 4. In sect. 5 we discuss the three 
dimensional reduction of the BS-equation and its application to the determination 
of the spectrum and the Regge trajectories. 

2 Green functions and Feynman-Schwinger rep- 
resentations 

The quark-antiquark and the single quark gauge invariant Green functions are de- 
fined as 

1 — — c 

G^\xi,X2,yi,y2) = -{0\Tij2ix2)U{x2,xi)i)i{xi)ip^{yi)U{yi,y2)i'2iy2)\0) = 

= -TTc{U{x2,xi)Si{xi,yi;A)U{yi,y2)S2{y2,X2;-A)) (3) 

(plus an annihilation term in the equal fiavor case) and 

G^\x -y) = (0|T?7(y,x)^(a;Mi/)|0) = ii:Tc{U{y,x)S{x,y-A)) , (4) 

where ijj'^ denotes the charge-conjugate fields, the tilde and Trc the transposition 
and the trace respectively over the color indices alone and U the path-ordered gauge 
string (Schwinger string) 



U{h, a) = P exp \ ig / dx'' A^{x) \ 



(5) 



Notice that the integration in (^ is along an arbitrary hne joining a to 6 (which 
usually we shall not specify explicitly), S, Si and 5*2 are the quark propagators in the 
external gauge field Af^ and the angle brackets denote average on the gauge variable 
alone (weighted in principle with the determinant Mf{A) resulting from the explicit 
integration of the fermionic fields). 

The propagator S is supposed to be defined by the equation (we shall suppress 
indices specifying the quarks, as a rule, when dealing with single quark quantities) 

{irD,-m)S{x,y;A)=6\x-y) (6) 

and the appropriate boundary conditions. This can be rewritten as JIl 

S{x, y- A) = {lYD, + m)/V{x, y- A), (7) 

in terms of a "second order" propagator defined in turn by the equation 

(D.D^ + m'~^g a^'^^F^,) A'^(x, y; A) = -6\x - y) , (8) 

with a'"' = ^[-f'',Y]- 

After replacing (|^) in (|^) and (|[), using an appropriate derivative it is possible 
to take the differential operator out of the angle brackets and write Q 

G^\xi,X2;yi,y2) = -{hidif, + mi){i'j!^d2u + m2)H'^\xi,X2;yi,y2) , (9) 

G^\x -y) = {i^'td, + m)H^\x - y) , (10) 

with 

H^\xi, X2-, 2/1, 2/2) = -^Trc(f/(a;2, Xi)A^(xi, y^- A)U{yi, 2/2)A^(x2, 2/2; -A)) , (11) 

H^\x -y) = zTrc(f/(2/, x)A'^(x, y- A)). (12) 

For the second order propagator we have the Feynman-Schwinger representation 

-^ fOQ rx rs ^ rs 

A''{x,y;A) = -- ds Vzexp[-i dT-{m'^ + z^)]S^F exp[ig drz^AJz)], 
2 JO Jy Jo 2 JO 

(13) 
with 

5S = Texp[-i/;...-^], (14) 



^ Given a functional $[7ab] of the curve ^ab with ends a and 6, let us assume that the variation 
of $ consequent to an infinitesimal modification of the curve 7 ^ 7 + ^7 can be expressed as the 
sum of various terms proportional respectively to 5a, to 5h and to the single elements 5S^'^{x) 
of the surface swept by the curve. Then, the derivatives daP, dbp and 5/5S''"{x) are defined by 
the equation 5$ = SaPdaP^ + SbPdbP^ + 5/ SSP'^{x)S^/SSP''{x). For a Schwinger string we 

have SU{b,a) ^ 5hPigAp{h)U{b,a) ~ 5aPU{b,a)igAp{a) + f jUSP''{z)P{-Fp„{z)U{b,a) and so 
BapU = -tgUApia), BbpU = igAp{b)U and jg^^U = V[-igFp,{z)U]. 
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T and P being the ordering prescriptions along the path acting on the spin and on 
the color matrices respectively and 5S^^ = dz^5z'^ —dz'^5z^ (the functional derivative 
being defined through an arbitrary deformation, z —^ z + 5z, of the line connecting 
a to 6, see footnote) [0]. Replacing ( |T3|) in ([TTD and (0) we obtain 



1 p^^ f^^ P"^^ /"^^ 

H^\xi,X2;yi,y2) = (t^)^ / dsi / ds2 / Vzi / Vz2 

I Jo Jo Jyi Jj/2 

exp{--y dTi{ml + zl)--J dT2{ml + zi)} 
S^^S',^]-{TTFexp{zgi dz^A,{z)}), (15) 



ilf^Xx-y) = -/ ds Vzexp{- - drirn^ + z^) 

2 Jo Jy 2 JO 

5o^(TrPexp{^^/ dz^A^{z)]). (16) 

Here, the loop Tqq occurring in the 4-points function is made by the quark world 
line 7i, the antiquark world line 72 followed in the reverse direction, and the two 
strings 0:1X2 and 2/21/1; the loop Vq occurring in the 2-points function is made simply 
by the quark trajectory 7 connecting y to x and the string yx. 

3 Wilson loop correlators 

We now apply Eq. (Q) to evaluate the Wilson correlators for Vqq and Vq. 
At the lowest order the perturbative term can be written for any loop F 

^(lniy[r])pert = -'^g''jdz^jdz'"D^,{z - z') (17) 

D^u{z — z') being the free gauge propagator. If we neglect the contribution coming 
from propagators connecting a point on a world-line to a point on a string or two 
point on the strings, we can write for Tqq 

4 „ /■*! /""a 

"2 



4 /-SI /•S2 

i(lnTy[rgg])pert = -9 dTi dT2D^^{zi~ Z2)z^zl 

o Jo -/ 

--/ E f ' ^^. /7 ^^;^.^(^. - 4)^'^r (18) 

and for F^ 

^(lniy[F,])pert = -^^' /'rfr r dT'D,,{z - z')z^z"' (19) 

O JO JO 



Let us further consider for the moment the case in which Tgg hes on a plane. Then 
^min coincides simply with the portion of plane delimited by Tgq and it can be 
written, in a four dimensional language, []1| 

Smin = / dn dT26{zio - Z2o)\zi - Z2\e{zio)e{z2o) 
Jo Jo 

/ dA|i;^o^20 - ('^ZiTi20 + (1 - A)z2Ti:io)^}'' - 
- Yl / dTj / dr'jSizjo - z'jo)\zj - z'j\e{zjo)e{z'jf^) 

7 = 1 ^ '^ 

r d\{z%z'^ - (AziTi^o + (1 - A)z;.Ti,o)'}' , (20) 

where we have used Zj for 2;j(tJ) and e(t) denotes the sign function. Eq. (pOD 
corresponds to span the surface by a straight line joining two points with the same 
time coordinate on the quark and the antiquark world lines respectively. The sign 
factors and the second term are necessary to reconstruct the surface as the algebraic 
sum of various pieces when the world-lines go backward in time. 

The "straight line equal time approximation" consists in assuming ( pOD even 
if Tgg does not stay on a plane. Notice that ( PU[ ) gives always 6*111111 correctly up 
to the order (z/io)^ in a semi-relativistic expansion. Furthermore (pOD is exact 
for particular geometries. This is the case e. g. for 71 and 72 making a regular 
double helix (with axis parallel to the time axis) corresponding to a pure rotational 
motion of the quark and the antiquark around a fixed point. Notice also that the 
approximation depends in general on the reference frame. Keeping in mind the helix 
example, however, we shall assume its validity in the center of mass frame. 

In a similar way, in the case of Tg we can set for x° = y^ we can set 

5*111111 = - dr dT'6{zo- ZQ)\z-z'\e{zo)e{zQ) dXizlz'^ - 
Jo Jo Jo ^ 

-(AzT4 + (l-A)z^io)'}'. (21) 

It is clear that, to include consistently the cases x^ 7^ X2, y\ 7^ ^25 ^° 7^ 2/°5 
the line integrals in (PDj ) and (^) should be extended in an obvious way to the 
Schwinger strings. However, we find convenient to consider two new functions that 
are obtained replacing (|l^) and ( pO|) in ( |T5l) and (|l^) and ( ^If ) in (|16]) so as they 



stand. We obtain in this way the following equations 

]^ ^CXD /-OO pX\ PX2 

H{xi,X2;yi,y2) = (7^)^ / dsi / ds2 / Vzi / Vz2 

I Jo Jo Jyx Jy2 

exp { - ^ $: fj dT,{m] + i|)}5o^^5o^^ exp [t Y^ £ dr^ JJ dr'^E{z, - z'^ i^, z'^) 

p dn P dT2E{z^ - Z2] ii, Z2)} (22) 



/o 



and 

][ roo rx Q rs 

H{x-y) = - ds Vzexp{-- dTim^ + z"^)} 

2 Jo Jy 2 JO 

S'^expiz r rEiz-z';z,z')} , (23) 

JO JO 

where we have set 

EiC;p,p') = ^pert(C;p,p') + ^conf(C;p,p') (24) 

with 



Epert = 47r|asD^.(C)pV 

i^conf = 5{Co)\chpoHp'o) Io' d\{pip's - [Ap[,PT + (1 - A)popyn^ 



(25) 



Obviously for arbitrary arguments the quantities H{xi,X2; yi, ^2) and H{x—y) as 
defined by (p^)-(p5D can differ very significantly from the original H^\xi,X2',yi,y2) 
and H^\x — y). However, as we mentioned, the two couples coincide in the limits 
X2 — >■ Xi, 1/2 — ^ Z/i, and y ^ X and they are completely equivalent for what concerns 
bound state problems and condensate determination. 

4 Bethe-Salpeter and Dyson-Schwinger equations 



From eq.s (22) and (^), by various manipulation and using an appropriate iterative 
procedure, a Bethe-Salpeter equation for the function H{xi, X2] yi, 1/2) and a Dyson- 
Schwinger equation for H{x — y) can be derived in the form |l| 

H{xi,X2;yi,y2) = Hi{xi - yi) H2{x2 - y2) ~ 

-i I d^C,id^^2d\d\Hi{xi - C.i) ^^2(3^2 - 6) 

x4fe(6, 6; ^1, V2) < (^2 Him, V2; yi, 1/2) , (26) 

H{x-y) = Ho{x-y)+i fd^^d^7]d^CdWHo{x-0 

x4,(e, e'; V, v'KHiv - r^')a'H{i' - y) , (27) 

where we have set a, 6 = 0, /ii/, with a'' = 1, and Hi and if 2 denote the quark and 
the antiquark H-propagators respectively. 

If we pass to the momentum representation, the corresponding homogeneous 
BS-equation becomes in a 4 x 4 matrix representation 

ffi(ip + k)a'1>p(u)a''H2{ - ip + *:) , (28) 



where $p(/c) denotes an appropriate wave function and the center of mass frame 
has to be understood; i.e. P = (m^, 0). 

Similarly, in terms of the irreducible self-energy, defined by H{k) = HQ{k) + 
iHo{k)r{k)H{k) , the DS-equation can be written also 

d^l f /, , k + l k + l^ _^,,„ ^ 



T{k) = j -^^h,[k-l-^^,^^)a'^H{l)ar (29) 

Notice that in principle ( PB| ) and ( P7| ) or ( P^D and (^Uj) are exact equations. 
However the kernels lab are generated in the form of an expansion in a^, and a. At 
the lowest order in both such constants, we have explicitly 

+^<y j £Ce~''^<\Q\e{p,)e{p',) f^ d\{plp'^ - [Ap^^pT + (1 - \)pov'tYV^ 






U,p.{Q;p,p') = -Am^a,p''{6lQ, - 6lQp)b^p{Q) + 



^ ICIv/p?-PT 



W-(Q;p,p') = 7r^as(5^Q. - 5^Q^)(5,"Q. - 5:Qp)D„^(Q) (30) 

where in the second and in the third equation Co = has to be understood. 
Setting 



I 



iH-\k) = J2Mk)hrik), (31) 



r=0 



with uq = 1, uji = 7°, ci;2 = —7 ■ k, Ws = 7^7 ■ k, k = -j^k and hQ^k), . . . h^^k) 
functions of /cq and |k|. Eq.(p^) can also be written 



RUk, I) = T\iab{k -l;^, ^)TT[u;nkK^s{k)a'] , (33) 

where the sign — applies to the s = case, the sign + to all the other cases. Notice 
that actually only Rqq, Ru, R12, R21, R22, R33 are different from zero. 

In connection with the problem of the light pseudo-scalar mesons, let us now 
consider eq.s (E3) and (p2D in the chiral limit mi = 1712 = 0. As it is apparent from 



(P^) chiral symmetry is broken in such limit, if hi and /i2 do not both vanish. On 
the other side, if this is the case and if in addition h^ ^ 0, it can be checked that 
eq. (0) is solved for P = by <l>o(A;) = Hi{k)[-f^h{k) +7 cdotkh2{k)]-f^H2{k), 
and a zero mass pseudo-scalar bound state exists, consistently with the Goldstone 
theorem. Notice that such a result is strictly related to the occurrence of the same 
kernel in the two equations and this in turn is due to the inclusion of the self-energy 
terms in (pO]), i. e. to the correct account of world lines that go backwards in time. 
Obviously the supposed behavior of the solution is superficially consistent with the 
form of (P2]); presently, however, we are not able to produce any proof. 



5 Three-dimensional reduction and spectrum 

By replacing Hi{k) and H2{k) in the BS equation by the corresponding free prop- 
agators ^2!^! and performing an instantaneous approximation on the kernels, one 

j 
can obtain a three dimensional reduction of the original equation in the form of the 

eigenvalue equation for the relativistic Hamiltonian |jl| 



H 



m^ 



k2 



m 



^ + k2 + K 



(34) 



with 



(k|^|k') 



2JwiW2w'iW 



4a« 



i,,/^ I3 7r2 



+ :^^' X k- {(Ti + (T2) + 



2Q2 



+ 7: CTi •cr2 





1 f 2 (Q-q)'\ 



- - 

4 V3 



g2o(«i ■ Q) -gio(«2 ■ Q) 
1 



CTl ■ 0-2 — 



2Q2 

(Q-o-i)(Q-o-2 



Q^ 



(27r)= 



rf're**^-V-^Hr,q,gio,g2o) , 



+ 



4Q2 



ai-Q)(«2-Q) 



(35) 



/°^*(r,q,gio,g2o) 



ar 



qio + <?2o 



?2oV?io - Qt + QwyQw - Qt 



,^10^20/ • IQtI . . IqtI 
-| — ; arcsm \- arcsm 



IQtI 



?io 



Q20 ~ cLt 



QW C[20 

r X q- cr2 - iq2o{r ■ 0L2 



'?20 



?io - Qt 



r X q- cTi + zgio(r • otx 

(36) 



61) k' and k denote the final and the initial center of mass momentum 



In eq.s 

of the quark; Wj = Jraj + k^, Wj = Jmf+k^, q = ^^ , Q = k' — k, Qjq 



Wj+w'. 



qt^ = {6 — f f'^^q'^ is the transverse momentum, a are the usual Dirac matrices 
7^7'^, and a'^ = i/2e''""^a"'"^ the 4x4 spin matrices. 
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Figure 1: Charmonium spectrum. 
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In spite of its complication the above expression has various significant hmit 
cases that corresponds to models successfully used in different areas. In the static 
limit it gives the local potential 



V = h err . 

3 r 



(37) 



In the heavy masses limit, by an — expansion and an appropriate Foldy-Wouthuysen 
transformation, it reduces to the potential discussed in ref . p| , If the spin dependent 
terms are neglected, V becomes identical (apart from a question of ordering) to the 
first order expansion of the potential corresponding to the relativistic fiux tube 
model 0, §. 

We have attempt to apply the actual V as given by (|3^) and (|36|) to the deter- 
mination of the spectrum. In the preliminary calculation we have performed up to 
now the spin orbit terms have been omitted, due to their complication; however, the 
hyperfine terms have been included. 

The numerical procedure we have followed consists in solving first the eigenvalue 
equation for the static potential ( ^7\j by the Rayleigh-Ritz method and then in 
evaluating the quantities {H) for the eigenf unctions obtained in the first step. We 
have adopted the following parameters: Og = 0.363, a = 0.175 GeV^, rUc = 1.405 
GeV, mb = 4.81 GeV, m, = 200 MeV, m„ = 10 MeV. The first four values have 
to be compared with those obtained from heavy quarkonium fits (e.g. 0, |^) and 
apart from the possibility of a small rearrangement are completely determined by 
these. The light quark masses have been fixed a priori on typical current values as 
reported by the Particle Data Group. 

By such parameters one succeeds to reproduce reasonably well the not only the 
bottonium and the charmonium spectrum, but also the Regge trajectories (with 
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7 M' 

ao(1450) 
Figure 2: Ground triplet uu Regge trajectories. Theoretical results (full line) com- 
pared with experimental data (circlet). Cross denote less established masses. 



Table 1: Theoretical results for uc, ub, sc, sb systems (MeV). Experimental data 
are enclosed in brackets. 



State 



uc 



ub 



sc 



sb 



IS 1973 (1973 ± 1) 5326 (5313 ± 2) 2080 (2076.4 ± 0.5) 5418 (5404.6 ± 2.5) 

2S 2600 (2623±?)^ 5906 (5897±?)^ 2713 6004 

IP 2442 (2438±?)^ 5777 (5825 ± Uy 2528 (2535.35 ± 0.34) 5848 (5853 ± 15) 

^Obtained from preliminary Delphi data m{D*') = 2637±8 MeV, m{B*') = 5906±14 

MeV subtracting 1/4 theoretical hyperfine splitting reported in table Q. 

^Estimated from m{D*) = 2459 ± 4 MeV, m{Di) = 2427 ± 5 MeV. 

'^From preliminary Delphi data R]. 



Table 2: Theoretical results for qq hyperfine splitting (MeV). Experimental data are 
enclosed in brackets. 



State 



uc 



ub 



cc 



bb 



sc 



sb 



IS 111 (141 ±1) 59 (46 ±3) 97 (117 ±2) 102 108(144) 60 (47 ± 4) 
2S 59 38 59 (92 ± 5) 42 62 40 
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correct slope and intercepts) for the ground triplet states of the uu, us, ss systems, 
and the known spin averaged states for the light-heavy systems. Notice that no ad 
hoc constant has been added to the potential and that in particular the heavy-light 
sector is completely parameter free. 

Our results have been reported in full in ref 0], where even some details on 
the the numerical difficulties are explained. Here as an example for the heavy-heavy 
systems we report in fig. |1] the charmonium spectrum and for the light-light systems 
we report in fig. ^ the p Regge trajectory. The results for heavy-light systems are 
reported in table |I| and compared with the experimental spin averaged masses by 
using the theoretical splitting of table ^ where the singlet states are not available. 
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